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Abstract. The newly proposed combinational approach, called the grand combinational model (GCM), as 
would be described in detail in the text, is still under our careful scrutiny. By applying it, we have attempted 
to analyze here the characteristics of both the transverse momentum(pT)-, and centrality-dependence of 
production of the main varieties of the secondaries measured in AuAu collisions at BNL-RHIC at both 
Y^snn = 130 GeV and ^/snn = 200 GeV by PHENIX Collaboration. Besides, with the help of it, we have 
also investigated the nature of centrality-dependence of the average transverse momenta of the various 
major categories of particles in AuAu collisions at RHIC. The model seems to survive quite smoothly the 
acid tests of the latest PHENIX data, as it accommodates data modestly well on these twin aspects. The 
study reveals a kind of universality of nature of the hadronic secondaries and also of the basic particle and 
nuclear interactions at high energies. However, in the end, we precisely point out both the strengths and 
limitations of the specific model under consideration here. 

PACS. 13.60.Hb Inclusive cross-section - 25.75.-q Relativistic Heavy Ion Collision. 



1 Introduction 

In the recent past the PHENIX Collaboration conducted 
exhaustive measurements of the particle-yields in AuAu 
collisions at RHIC-BNL at very high ranges of transverse 
momentum values and also at various centralities of the 
collisions. The data-sets on both identified centrality- de- 
pendent charged hadron spectra vs. transverse momenta 
and also on the average transverse momentum- values (< 
Pt >) of the various secondaries vs. the centralities of 
the collisions'^] were indicated by the graphical plots 
shown in various figures presented here. These measure- 
ments were mainly aimed at checking some predictions 
from the Standard Model(SM) on (i) the suppression of 
particle yields at large transverse momenta and (ii) the 
dependence of the degree of suppression on the centrality 
of the collision. In the framework of the standard model, 
all the secondary particles are produced by jets. These 
jets are likely to suffer significant energy loss via 'gluon' 
radiation, while passing through the hot dense medium 
arising out of the anticipated QGP formation. So, accord- 
ing to this view, the particle yield at large pt should be 
suppressed to a considerable degree. Secondly, since the 
amount of energy loss is a function of the density and the 
path length through such a projected QGP media, the 
suppression effect is predicted to be dependent also on 
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the centrality of the collision^]. All the theoretical pos- 
tulates and predictions serve as the background for such 
intense attempts at measurements. 

But in our approach we would maintain a degree of 
neutrality to such standard theoretical views and would 
concern ourselves mainly with the results of measurements 
with the ulterior motive of checking and testing an alter- 
native approach, called here the Combinational Approach 
or the grand combination of models(GCM) which would 
be utilized here to explain the sets of data on: (i) identified 
charged hadron spectra vs. pi-values and (ii) the values 
of average transverse momentum, denoted by < >, 
vs. the different centrality-ranges of the AuAu reaction in 
particular. 

The work is organized as follows: In Section 2 we offer 
an outline of the GCM to be made use of in the present 
study. The Section 3 contains and depicts the concrete 
results obtained on the basis of this model in the tabular 
forms and the graphical plots. The last section presents, 
as usual, the summary and final comments. 



2 The Model: A Sketch 

Following the suggestion of Faesslerpj] and the work of 
PeitzmannJS] and also of Schmidt and Schukraft 6 , we 
propose here a generalized empirical relationship between 
the inclusive cross-sections for any variety of the secon- 
daries (Q), such as hadrons, pions, kaons or proton/ antiprotpons, 
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produced in nucleon(N)-nucleon(N) collision and that for 
nucleus (A)-nucleus(B) collision as given below: 



QX) 



( j4B )9(Pt) rpp _> QJQ 

dp A 



The main bulk of work, thus, converges to the mak- 
ing of an appropriate choice of form for g(px)- And the 
necessary choices are to be made on the basis of certain 
premises and physical considerations which do not violate 
the canons of high energy particle interactions. 

The expression for inclusive cross-section of Q in proton- 
proton scattering at high energies occurring in Eqn.(2) 
could be chosen in the form suggested first by HagedornJJ]: 



QX) = Ci ( 1 + 



PQ ' 



(3) 



where C\ is the normalization constant, and p Q , n are 
interaction-dependent chosen phenomenological parame- 
ters for which the values are to be obtained by the method 
of fitting the spectra in PP interaction. 

The final working formula for the nucleus-nucleus col- 
lisions is now being proposed here in the form given below: 



precise determination of e is not possible for the following 
understated reasons: 

(i) To make our point let us recast the expression for 
(4) in the form given below: 



QX) ~ (AB)+w> p ^ £— r (PP -» QX) , 
dp 6 

(1) 

where <fr(y, pt) could be expressed in the factorization 
form, <p(y, px) = f(y) g{pr)', and the product, AB on the 
right hand side of the above equation is the product of 
mass numbers of the two nuclei participating in the colli- 
sions at high energies, of which one will be the projectile 
and the other one the target. 

While investigating a specific nature of dependence of 
the two variables (y and pt), either of them is assumed 
to remain averaged or with definite values. Speaking in 
clearer terms, if and when px-dependence is studied by 
experimental group, the rapidity factor is integrated over 
certain limits and is absorbed in the normalization factor. 
So, the formula turns into 



E 



^{AB -> QX) oc (4B)('+»pt-* t ) e^(PP -> QX) 
cx (AB)( e + a P T - ^t) (1 + £2L)-n 

(4) 

with <?(pt) = (e + apT — (3p\), where this suggestion 
of quadratic parametrization for q(pt) is exclusively made 
by us and is called hereafter De-Bhattacharyya parametriza- 
tion(DBP). In the above expression e, a and (3 are con- 
stants for a specific pair of projectile and target. 

Earlier experimental works 8 9, 10 showed that g(pr) 
is less than unity in the px-domain, px < 1-5 GeV/c. 
Besides, it was also observed that the parameter e, which 
gives the value of g{pr) at px = 0, is also less than one 
and this value differs from collision to collision. The other 
two parameters a and j3 essentially determine the nature 
of curvature of g(j?x)- However, in the present context, 



QX) sa C 2 (AB) e {AB) ( > ap T- (lp ^(l + ?£.) 



Po 



where C2 is the normalization term which has a depen- 
dence either on the rapidity or on the rapidity density 
of the Q and which also absorbs the previous constant 
term,Ci as well. 

Quite obviously, we have adopted here the method of 
fitting. Now, in Eqn.(5) one finds that there are two con- 
stant terms C2 and e which are neither the coefficients 
nor the exponent terms of any function of the variable, 
Pt. And as e is a constant for a specific collision at a spe- 
cific energy, the product of the two terms C2 and (AB) e 
appears as just a new constant. And, it will just not be 
possible to obtain fit-values simultaneously for two con- 
stants of the above types by the method of fitting. 

(ii) From Eqn.(2) the nature of g(px) can easily be 
determined by calculating the ratio of the logarithm of 
the ratios of nuclear-to-PP collision and the logarithm 
of the product AB. Thus, one can measure e from the 
intercept of g(pr) along y-axis as soon as one gets the 
values of -E/^-f for any specific secondary production in 
both AB collision and PP collision at the same cm. en- 
ergy. In the present study we have tried to consider the 
AuAu collision system in various centrality bins at two 
different c.m.energis. In order to do so, we have to con- 
sider the data on normalized versions of E^^ for any sec- 
ondary particle produced in this collision system for which 
no clear £J^-f-data is available to us. Furthermore, from 
these normalized versions we can/could not extract the ap- 
propriate values of E^- as the normalization terms, total 
inclusive cross-sections(c7i n ) etc., for this collision system 
at all centrality-bins cannot always be readily obtained. 
Besides, it will also not be possible to get readily the data 
on inclusive spectra for PP collisions at all cm. energies. 

In order to sidetrack these difficulties and also to build 
up an escape-route, we have concentrated here almost 
wholly to the values of a and (3 for various collision sys- 
tems and the resultant effects of C2 and e have been ab- 
sorbed into a single constant term C3. Hence, the final 
expression becomes 



E 



QX) » C 3 (AB) 



(ap T -/3p T ) 



( 1 



Pt_ 

Po 



r 



(6) 

with C 3 = C 2 (AB) e . 

The exponent factor term cvpx — @p\ obviously repre- 
sents here [g(pr) — e] instead of <?(px) alone The expres- 
sion^) given above is the physical embodiment of what we 
have termed to be the grand combination of models(GCM) 
that has been utilized here. The results of PP scattering 
are obtained in the above on the basis of eqn.(3) provided 
by Hagedorn's model(HM); and the route for converting 
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the results of NN to NA or AB collisions is built up by the 
Peitzmann's approach(PA) represented by expression(2). 
The further input is the De-Bhattacharyya parametriza- 
tion for the nature of the exponent. Thus, the GCM is 
the combination of HM, PA and the DBP, all of which are 
used here. 

And the choice of this form of parametrization for the 
power of the exponent in eqn.(4) is not altogether a co- 
incidence. In dealing with the EMC effect in the lepton- 
nucleus collisions, one of the authors here(SB').|ll| made 
use of a polynomial form of A-dependence with a variable 
x which is a variant of xf (the Feynman Scaling vari- 
able). This gives us a clue to make a similar choice for 
both g(pr) and f(y) variable(s) in each case separately. 
In the recent times, De-Bhattacharyya parametrization is 
being extensively applied to interpret the measured data 
on the various aspects [T2*lll3III4III5| of the particle- nucleus 
and nucleus-nucleus interactions at high energies. In the 
recent past Hwa et. al.^H] also made use of this sort of 
relationship in a somewhat different context. The under- 
lying physics implications of this parametrization stem 
mainly from the cxprcssion(4) which could be identified as 
a clear mechanism for switch-over of the results obtained 
for nucleon-nucleon(PP) collision to those for nucleus- 
nucleus interactions at high energies in a direct and straight- 
forward manner. The polynomial exponent of the product 
term on AB takes care of the totality of the nuclear effects. 

For the sake of clarity and confirmation, let us fur- 
ther emphasize a point here very categorically. It is to be 
noted that this model(GCM) containing all the Eqns.(4), 
(5) and (6) was described in some detail earlier and was 
made use of in analyzing extensive sets of data in the pre- 
vious publications 1 1 211 1 311 1 5) by the same authors. And in 
verifying the validity of this model further, the purpose 
here is to apply the same model to some other problem- 
atical aspects of data which we would dwell upon in the 
subsequent sections. Before taking them up, let us state a 
point. In some previous works [T3lll4lll5j we tried to pro- 
vide some sort of physical interpretations for some of the 
parameters used in the present work. But, those explana- 
tions were only of suggestive nature. Besides, obviously, 
they are not complete and sufficient, for which we have 
chosen not to reiterate them here once more. 



3 Presentation of Results 



Obviously the GCM is the model of our choice here. The 
inclusive spectra for production of the main varieties of 
various secondaries produced in AuAu collisions at RHIC 
at both v / snn = 130 GeV and ^/snn = 200 GeV and also 
at various centrality values have been worked out here 
phenomenologically and shown in the several diagrams. 
The values of po and n, occurring in eqn.(3), which are es- 
sentially the contribution of PP collisions to the nucleus- 
nucleus collisions for the same secondaries produced at the 
same cm. energies per nucleon, have been introduced by 
the following relationships: 



(7) 



(8) 



These are products of just empirical analyzes made 
earlier and reported in some of our previous works [T2l 
II3| . The actually used values of the arbitrary parameters, 
a,b,a and b for various secondary particles are given in 
Table- 1. 

The obtained values of the average yields of hadrons at 
various centralities have been depicted in Fig.l. The left 
panel is for a/snn = 130 GeV and the right panel is for 
v / snn = 200 GeV. The Fig. 2 describes pion production 
on a charge-neutral and average basis at both the ener- 
gies. Similar is case with kaons in Fig. 3. The diagrams in 
Fig. 4 reproduce data on proton production in AuAu colli- 
sions at the two RHIC energies. The solid curves in Fig. 5 
display the theoretical yields of antiprotons in AuAu re- 
action at various values of the centrality of the collision 
against measured data. The values of the parameters a 
and P have been initially chosen with the singular moti- 
vation of obtaining satisfactory fits to the data, though 
finally even this arbitrariness has led to some revelation 
of the specific nature of a and (3 as are shown in the plots 
of Fig. 6 for the various secondaries in several panels. 

The values of a and (3 to be used in obtaining our 
model-based results are shown in the different tables(Table 
2 - Table II). The extreme left columns in all of them con- 
tain information about the centrality of the reaction and 
the extreme right ones offer the x 2 / ncu " values. Even for 
the cases of proton and antiproton production, wherein 
the data suffer a high degree of uncertainty, the x 2 /ndf 
values are modestly satisfactory. The systematic trends 
of the used values of a and (5 depict a harmony of their 
nature which have been hinted by Table 12 and Table 13 
and represented by the sets of diagrams in Fig. 6. The solid 
lines in Fig. 6 provide the phcnomcnological fits which can 
be expressed by a common relationship of the form given 
below, 



a(A par t), /3(A p 

art / 



I? 



S In (A part ) 



(9) 



The different values of R and S for various secondaries 
are given in Table 12 and Table 13. 

A comment is in order here in a preemptive manner. 
The values of a and /3 shown in our previous work |12j even 
on AuAu collision could be and are little different from 
what are depicted here for the two reasons: (i) The pt- 
range of the detected secondaries in the previous work was 
limited mostly in the region from 0.8 GeV/c to 3 GeV/c, 
whereas in the present case both the full low pt and a 
larger domain of high p^ range for the secondaries (mainly 
charged hadrons) has been covered, (ii) Secondly, in the 
former study [T2*] the minimum bias event was studied in 
the main. On the contrary, the present study has very 
much been centrality-specific and the data for the various 
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centrality-values of the AuAu collision have been served 
within a phenomenological framework. 

The average transverse momenta values for the differ- 
ent categories of particles have, however, been worked out 
on the basis of the following expression: 



<Pt > 



ft 



(IN 



p T 



dpL 



(10) 



The values of a and (3 to be introduced for < pt >- 
values are used in both energy-specific and particle-specific 
manner with the help of eqn.(9). The GCM-based results 
on < pt > values are plotted in Fig. 7. The different 
centrality values, the particle-species and the interaction 
energy- values are separately mentioned in each of the di- 
agrams. 



4 Concluding Remarks 

The chosen model appears to present essentially a uni- 
versal approach in the sense that (i) it provides a uni- 
fied description of data on particle production in nuclear 
collisions in terms of the basic PP interaction; (ii) the 
method could be applied in an integrated and uniform 
way without introduction of any artificial divide between 
the so-called 'soft'(low-pr) and 'hard'(large-pT) interac- 
tions; (iii) the general approach remains valid, irrespec- 
tive of whether the collisions are central or peripheral; 
(iv) it has no model- or mechanism-specific physical pic- 
ture as the input and as the constraint as well; (v) the 
values of a and (3 which are the only arbitrary parame- 
ters need to be assumed and they demand tuning and ad- 
justment on a case-to-case basis in an interaction-specific, 
secondary-specific and centrality-specific manner. So, the 
model seems to provide a universal, useful and economi- 
cal description of a large body of data on the high energy 
heavy ion collisions. The model is useful, as it is seen to 
give a fair account of the vast amount of data; and it is eco- 
nomical, because there are only two arbitrary parameters 
alongwith one normalization term for the general studies 
on heavy ion reactions at high energies. 

The agreements between the measured and/or extracted 
data and the phenomenological outputs are quite satisfac- 
tory on an overall basis of the px-spectra and the < pt > 
vs. centrality diagrams. The only exception is the case of 
the average transverse momenta values of kaons in AuAu 
reaction at ^/snn = 200 GeV. Though we cannot readily 
ascribe any reason for such departure, there is a general 
observation that the measurements related to any variety 
of the strange particles suffer, in general, a higher degree 
of uncertainty. Besides, we also fail to explain here how 
and why this phenomenological approach works function- 
ally so well. No clue to any concrete physical reason arising 
out of some underlying dynamics of particle and nuclear 
interactions could be immediately ascertained. The har- 
mony revealed in the natures of a and (3 versus various 
^part-values reflecting various centrality of the reactions 
is certainly an interesting observation from the present ap- 
proach. The closeness of the values of S in the Table- 12 



and Table- 13 at two different neighbourly energies appears 
to indicate the fact that, beyond a certain value of iV par t, 
the enhancement of the centrality of the collision with the 
increase in the number of wounded nucleons, i.e. N paxt 
does not necessarily and appreciably raise the values of 
< Pt >s for any variety of the secondaries. This brings 
out strong hints to what is called parton saturation at 
and after a definite value of A^ part . In fine, in so far as 
actual performance is concerned, the model has a modest 
degree of success. However, one major drawback in apply- 
ing this approach is it's over-reliance on the availability of 
the measured and dependable data-sets on the specific va- 
riety of the secondary in PP interaction at some definite 
energies and at certain reasonable intervals in the energy- 
values in order to construct the cncrgy-dcpcndcncc profile 
for some parameters to be used in the model. Secondly, 
the final working formula for studying the properties of 
nuclear collisions in the present work does neither contain 
directly, nor exhibit any of the technicalities of the nuclear 
geometry, e.g., the impact parameter (denoted generally by 
b) or of the space-time evolution scenarios of the nuclear 
collisions. The entirety of the nuclear effects is taken care 
of by the simple product term (AB)^ V,PT \ This simplic- 
ity of form could very well be viewed in a positive way in 
favour of the model. 



The authors would like to express their thankful gratitude to 
the anonymous referee for his/her critical remarks and valu- 
able suggestions for improvement of an earlier draft of the 
manuscript. 
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Table 1. Values of a,b,a' and b' used in eqn.(7) and in eqn.(8) to obtain po and n for various secondaries. 



Secondary Type 


a 


b 


a 


b 




1.5 


79.4 


6.5 


127 


K 


1.6 


103 


3.6 


161 


P 


7 


602 


5 


644 


P 


7 


478 


13 


527 



Table 2. Parameter values for charged hadrons(averaged) production in AuAu collision at ^/snn =130 GeV 



Centrality 


c 3 


a(c/GeV) 


/3(c/GeV) 2 


X7ndf 


0- 


-15% 


1708 ± 42 


0.10 ±0.01 


0.019 ±0.002 


1.127 


5- 


- 10% 


1487 ± 39 


0.09 ±0.01 


0.017 ±0.002 


1.281 


10 


-20% 


1161 ±31 


0.088 ± 0.004 


0.016 ±0.001 


1.401 


20 


-30% 


831 ± 21 


0.077 ±0.003 


0.012 ±0.001 


1.425 


30 


-40% 


561 ± 15 


0.073 ± 0.003 


0.012 ±0.002 


1.590 


40 


-60% 


262 ±5 


0.070 ±0.002 


0.011 ±0.001 


0.714 


60 


-80% 


69 ±2 


0.063 ±0.003 


0.009 ±0.001 


0.214 



Table 3. Parameter values for production of charged hadrons(averaged) in AuAu collision at ^/snn =200 GeV 



Centrality 


c 3 


a(c/GeV) 




xVndf 


0- 


-5% 


3015 ± 189 


0.072 ± 0.006 


0.012 ±0.002 


1.281 


10 


-15% 


1962 ± 120 


0.0741 ± 0.006 


0.012 ±0.002 


1.191 


20 


-30% 


1230 ± 111 


0.064 ± 0.009 


0.009 ± 0.002 


1.847 


40 


-50% 


547 ± 32 


0.052 ±0.004 


0.006 ±0.001 


0.943 


60 


-70% 


171 ±8 


0.041 ± 0.004 


0.005 ± 0.001 


0.403 


80 


-91% 


262 ±5 


0.070 ± 0.002 


0.011 ±0.001 


0.521 



Table 4. Parameter values for charged pions(averaged) produced in AuAu collision at ^/snn =130 GeV 



Centrality 


C 3 


a(c/GeV) 


/3(c/GeVr 


X>df 


0-5% 


1243 ± 113 


0.07 ±0.02 


0.022 ±0.010 


0.948 


5 - 15% 


922 ± 44 


0.07 ±0.01 


0.017 ±0.005 


0.399 


15 - 30% 


602 ± 47 


0.066 ±0.016 


0.017 ±0.008 


1.077 


30 - 60% 


230 ± 24 


0.062 ±0.022 


0.015 ±0.010 


1.774 


60 - 90% 


33 ± 10 


0.046 ± 0.010 


0.016 ±0.006 


1.859 



Table 5. Parameter values for charged pions(averaged) produced in AuAu collision at ^/snn = 200 GeV 



Centrality 


c 3 


a(c/GeV) 


[3(c/GeVy 


X Vndf 


0- 


-5% 


1375 ± 44 


0.068 ± 0.007 


0.024 ±0.003 


0.896 


10 


-15% 


1170 ± 26 


0.048 ± 0.005 


0.015 ±0.002 


0.850 


20 


-30% 


702 ± 24 


0.055 ±0.007 


0.017 ±0.003 


0.566 


30 


-40% 


502 ± 15 


0.039 ±0.006 


0.010 ±0.003 


0.817 


40 


-50% 


309 ± 14 


0.038 ±0.008 


0.012 ±0.003 


0.805 


50 


-60% 


181 ±4 


0.030 ± 0.004 


0.010 ±0.002 


1.159 


60 


-70% 


96 ±2 


0.022 ± 0.004 


0.010 ±0.003 


0.881 


70 


-80% 


44 ±1 


0.022 ± 0.006 


0.018 ±0.005 


1.623 


80 


-91% 


23 ± 1 


0.017 ±0.004 


0.010 ± 0.002 


1.753 



Table 6. Parameter values for charged kaons( averaged) produced in AuAu collision at ^/snn = 130 GeV 



Centrality 


c 3 


a(c/GeV) 




X 7ndf 


0-5% 


175 ± 38 


0.049 ±0.009 


0.018 ±0.007 


0.276 


5 - 15% 


135 ±6 


0.052 ± 0.009 


0.017 ±0.007 


0.489 


15 - 30% 


79 ±3 


0.047 ± 0.008 


0.016 ±0.008 


0.440 


30 - 60% 


32 ±4 


0.029 ± 0.006 


0.014 ±0.005 


0.150 


60 - 90% 


4.2 ±0.5 


0.019 ±0.008 


0.012 ±0.001 


0.322 
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Table 7. Parameter values for charged kaons(averaged) produced in AuAu collision at ^/ snn = 200 GeV 



Centrality 


Cz 


a(c/GeV) 


P(c/G<Nf 


X Vndf 


0- 


-5% 


150 ± 10 


0.081 ±0.010 


0.034 ±0.006 


0.622 


10 


-15% 


115 ±7 


0.080 ±0.011 


0.033 ± 0.006 


1.163 


20 


-30% 


70 ±9 


0.076 ± 0.009 


0.031 ± 0.005 


1.067 


30 


-40% 


45 ±8 


0.070 ± 0.009 


0.026 ± 0.006 


0.659 


40 


-50% 


30 ±7 


0.060 ± 0.008 


0.026 ± 0.005 


0.806 


50 


-60% 


15 ±5 


0.054 ± 0.007 


0.024 ± 0.005 


1.441 


60 


-70% 


7± 1 


0.043 ± 0.007 


0.022 ± 0.004 


0.996 


70 


-80% 


3.0 ±0.4 


0.034 ± 0.005 


0.021 ± 0.004 


1.762 


80 


-91% 


1.4 ±0.1 


0.025 ± 0.004 


0.020 ± 0.003 


2.015 



Table 8. Parameter values for production of secondary protons in AuAu collision at ^/snn = 130 GeV 



Centrality 


c. 


a(c/GeV) 


P(c/GeVy 


X 7ndf 


0-5% 


23 ±5 


0.21 ±0.03 


0.045 ±0.011 


1.028 


5 - 15% 


16 ±3 


0.22 ±0.02 


0.053 ± 0.008 


0.895 


15 - 30% 


12 ±2 


0.19 ±0.02 


0.042 ± 0.008 


0.936 


30 - 60% 


6.0 ±0.3 


0.15 ±0.01 


0.031 ± 0.003 


0.953 


60 - 90% 


1.0 ±0.1 


0.011 ±0.02 


0.023 ±0.006 


2.023 



Table 9. Parameter values for production of secondary protons in AuAu collision at ^/snn = 200 GeV 



Centrality 


c 3 


a(c/GeV) 


/3(c/GeVf 


X Vndf 


0- 


-5% 


19 ±3 


0.19 ±0.02 


0.037 ±0.004 


1.526 


10 


-15% 


17±4 


0.18 ±0.02 


0.034 ± 0.005 


0.927 


20 


-30% 


12 ±2 


0.16 ±0.02 


0.030 ± 0.004 


1.085 


30 


-40% 


8± 1 


0.15 ±0.01 


0.027 ± 0.003 


1.001 


40 


-50% 


6.5 ±0.7 


0.13 ±0.02 


0.025 ± 0.005 


1.703 


50 


-60% 


4.0 ±0.5 


0.12 ±0.01 


0.023 ± 0.003 


1.494 


60 


-70% 


1.9 ±0.3 


0.11 ±0.01 


0.020 ± 0.004 


1.896 


70 


-80% 


0.8 ±0.1 


0.10 ±0.01 


0.017 ±0.003 


2.116 


80 


-91% 


0.50 ±0.05 


0.08 ±0.01 


0.015 ±0.002 


1.839 



Table 10. Parameter values for production of secondary antiprotons in AuAu collision at snn =130 GeV 



Centrality 


c 3 


a(c/GeV) 


/3(c/GeV) 2 


XVndf 


0-5% 


13 ±2 


0.27 ± 0.02 


0.050 ±0.008 


1.336 


5 - 15% 


8±2 


0.28 ±0.03 


0.052 ±0.010 


1.652 


15 - 30% 


7.3 ±0.5 


0.24 ±0.01 


0.049 ± 0.005 


1.215 


30 - 60% 


3.7 ±0.5 


0.22 ±0.01 


0.046 ± 0.007 


1.658 


60 - 90% 


0.5 ±0.1 


0.18 ±0.02 


0.028 ±0.011 


1.958 



Table 11. Parameter values for production of secondary antiprotons in AuAu collision at ^/snn = 200 GeV 



Centrality 


c 3 


a(c/GeV) 


/3(c/GeVr 


XVndf 


0- 


-5% 


8.6 ±0.7 


0.27 ±0.02 


0.048 ± 0.003 


1.173 


10 


-15% 


8.5 ±0.8 


0.26 ±0.01 


0.046 ± 0.002 


1.102 


20 


-30% 


5.8 ±0.7 


0.25 ±0.01 


0.043 ± 0.003 


1.436 


30 


-40% 


4.1 ±0.6 


0.23 ±0.01 


0.042 ± 0.004 


1.269 


40 


-50% 


3.0 ±0.3 


0.22 ±0.01 


0.040 ± 0.003 


1.126 


50 


-60% 


2.2 ±0.3 


0.20 ±0.01 


0.039 ± 0.004 


1.059 


60 


-70% 


1.5 ±0.1 


0.17 ±0.02 


0.034 ± 0.003 


1.481 


70 


-80% 


0.6 ±0.1 


0.15 ±0.01 


0.033 ± 0.005 


1.470 


80 


-91% 


0.32 ±0.05 


0.14 ±0.01 


0.029 ± 0.006 


1.994 
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Table 12. Values of different parameters to obtain the 'Fit' for a on the basis of eqn.(9). 



Secondary type 


Collision energy 


/t(c/GeV) 


<b(c/Gev) 


2 


130 GeV 


0.036 


0.0085 




200 GeV 


0.018 


0.0085 


2 


130 GeV 


0.023 


0.0085 




200 GeV 


0.0003 


0.0087 


K ++K~ 
2 


130 GeV 


-0.011 


0.011 




200 GeV 


0.001 


0.015 


P 


130 GeV 


0.052 


0.026 




200 GeV 


0.034 


0.027 




130 GeV 


0.11 


0.025 


P 


200 GeV 


0.010 


0.026 



Table 13. Values of different parameters to obtain the 'Fit' for (3 on the basis of eqn.(9). 



Secondary type 


Collision energy 


R{c/GeVf 


S(c/GeV) x 


h++h- 

2 


130 GeV 


0.0035 


0.0021 




200 GeV 


0.0001 


0.0019 


2 


130 GeV 


0.008 


0.0019 




200 GeV 


0.0055 


0.0019 


K++K~ 
2 


130 GeV 


0.006 


0.0019 




200 GeV 


0.016 


0.0020 


P 


130 GeV 


0.012 


0.005 




200 GeV 


0.0039 


0.0048 




130 GeV 


0.022 


0.0051 


P 


200 GeV 


0.020 


0.0048 
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Fig. 1. Plots of invariant spectra of secondary charged hadrons produced in AuAu collisions at two different RHIC-energies 
for various centrality-bins. The experimental data points are taken from Ref.[I7 for ^/snn = 130 GeV and from Ref.0 for 
^/snn = 200 GeV. The solid curves provide the GCM-based fits. 
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Fig. 2. Nature of invariant spectra of secondary charged pions produced in AuAu collisions at two different RHIC-energies at 
different centralities as a function of px- The experimental data points at y / %N = 130 GeV are taken from Ref. 1]|] while those 
for ^/snn = 200 GeV are from Ref.0. The solid curves are drawn on the basis of eqn.(6). 
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Fig. 3. Transverse momentum spectra for production of secondary charged kaons in AuAu collisions at two different RH1C- 
energies at different centralities. The various experimental data points at ^/snn = 130 GeV are taken from Ref. 18. while those 
for -^snn = 200 GeV are from Ref.JT]. The present model(GCM)-based fits are depicted by the solid curvilinear lines. 



Au+Au 



130 GeV 




Au+Au, Vs NN = 200 GeV 



B (0-5%)x10, 
~ (5-15%)x10 2 , 
O (15-30%)x10 3 
(30-60%)x10, 
(60-92%)x10 7 
GCM 




Q (0-5%)x10 u 

• (10-15%)x10, 
O (20-30%)x10, 
m (30-40%)x10 3 
^ (40-50%)x10* 

* (50-60%)x10^ 
<V> (60-70%)x10^ 
w (70-80%)x10 8 „ 
^ (80-91 %)x10 10 
- GCM 



p T (GeV/c) 



(a) 



(b) 



Fig. 4. Plots of invariant spectra for production of secondary protons in AuAu collisions at two different RHIC-energies 
at different centralities. The various experimental data points at v / snn = 130 GeV are taken from Ref. 18 while those for 
y / %N = 200 GeV are from Ref.pQ. The GCM-based fits are depicted by the solid curves. 




Fig. 5. Nature of invariant spectra of secondary antiprotons produced in AuAu collisions at two different RHIC-energies at 
different centralities as a function of Pt- The experimental data points at ^/snn = 130 GeV are taken from Ref.|T2J while those 
for ^/snn = 200 GeV are from Ref.0- The solid curves are drawn on the basis of eqn.(6). 
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Fig. 6. Plots of a[in (c/GeV)] and [3[in (c/GeV) 2 ] as a function of number of participant nucleons, Af part . The different data-type 
points for various secondaries are taken from Table 2 - Table 11. Various solid curves are drawn on the basis of eqn.(9) and 
denoted as 'Fit' in the figures. 
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Fig. 7. Nature of average transverse momenta(< pt >) for various secondaries produced in AuAu collisions at two different 
energies as a function of 7V par t. The data-type points are the extracted-results obtained by RHIC experimental groupspp. The 
GCM-based results are shown by various dashed and dotted curves. 



